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ABSTRACT 

In this paper we report on the bulk features of the hole carved by the com- 
panion star in the material ejected during a Type la supernova explosion. In 
particular we are interested in the long term evolution of the hole as well as in 
its fingerprint in the geometry of the supernova remnant after several centuries 
of evolution, which is a hot topic in current Type la supernovae studies. We 
use an axisymmetric smoothed particle hydrodynamics code to characterize the 
geometric properties of the supernova remnant resulting from the interaction of 
this ejected material with the ambient medium. Our aim is to use supernova 
remnant observations to constrain the single degenerate scenario for Type la su- 
pernova progenitors. Our simulations show that the hole will remain open during 
centuries, although its partial or total closure at later times due to hydrodynamic 
instabilities is not excluded. Close to the edge of the hole, the Ray leigh- Taylor 
instability grows faster, leading to plumes that approach the edge of the forward 
shock. We also discuss other geometrical properties of the simulations, hke the 
evolution of the contact discontinuity. 

Subject headings: supernova: general, supernova remnants 
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Introduction 



The precise nature of tlie progenitors of Type la supernovae (SN la) is one of the most 
important unsolved problems in modern astrophysics. Because of the key role played by SN 
la in the chemical and dynamical evolution of galaxies and in cosmology the quest for these 
elusive progenitors has become a priority and a challenge to astronomers. Several lines of 
observational evidence point to the thermonuclear explosion of a white dwarf as the most 
probable progenitor because of the conspicuous absence of hydrogen emission lines in SN la 
spectra and the occurrence of these explosions in all galaxy types. Theoretical considerations 
also favour white dwarfs as SN la progenitors, because triggering a thermonuclear explosion 
in a degenerate object is not difficult, provided it has enough nuclear fuel to be burnt. 
However, the puzzle starts just beyond this point because there are many different ways to 
explode a white dwarf, each of them involving a different astrophysical scenario. Branch ^ 



Khokhlov ( 1995 ) , Hillebrandt & Niemeyer ( 2000 ) . 



One of these scenarios is the so-called single degenerate (SD) scenario, in which the 
explosion takes place in a compact binary system composed by a massive white dwarf and 
a non-degenerate secondary star, either evolved or unevolved. If the average mass transfer 
rate from the secondary is around ^ 10""^ M0.yr~^, surface ignition is avoided and the white 
dwarf manages to gradually increase its mass, eventually approaching the Chandrasekhar 



mass limit Hachisu, Kato & Nomoto (1996). At this point a carbon deflagration ensues close 



to the center of the white dwarf. Conductive and convective heat transport mechanisms 
spread the combustion to the whole mass of the star in a time scale of the order of a second, 
triggering a thermonuclear explosion. HydrodjTiamic models of this kind of explosion 
suggest that the observed nucleosynthetic yields and kinetic energies are better explained if 



the deflagration turns into a detonation (i.e. supersonic burning) at some point, Hoeflich & 



Khokhlov (1996), Gamezo, Khokhlov & Oran (2005), but a consensus theoretical model 
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for Type la supernovae does not exist yet. 



One way to constrain the identity of SN la progenitors is to assume that the SD 
scenario is generally valid, and then explore the observational consequences. Some of the 
implications of the SD hypothesis have to do with the effects that the presence of the 
nearby secondary star has on the bulk properties of the supernova ejecta immediately 
after the explosion, once the homologous stage has been reached. In this context, the 
presence of a companion star could modify the picture of the explosion in several ways, a) 
if sufficient material is stripped from the outer layers of the secondary during the collision 
with the ejecta, then this ablated, hydrogen rich, material should be seen in the SN spectra. 



Marietta, Burrows & Fryxell (2000). b) The large kick given to the secondary by the SN 
blast wave would dramatically change its internal structure and deposit a large amount of 
linear momentum. In this case, the observation of a peculiar star with large proper motion 
close to the expansion center of historical SN la remnants would lend clear support to the 
SD scenario, ( Canal, Mendez fc Ruiz-Lapuente||2001 ). c) The presence of the secondary star 
can also break the symmetry of the SN ejecta, introducing systematic effects that could be 



detected in spectral and spectropolarimetric SN la observations, Kasen et al. (2004). 



The first item above has been addressed in many works using both analytical 



calculations. Wheeler, Lecar & McKee (1975) and hydrodynamic models in two (Marietta 



Burrows & Fryxell 2000) and three dimensions (Serichol 2005 Pakmor et al. 2008). All 



these calculations agree that, depending on the nature of the companion, the stripped 
amount of hydrogen should range from 0.01-0.10 M0 for main-sequence (MS) stars to 
~ 0.5 M0 for red giants (RG), and therefore it should be detectable in SN la spectra . 
Taken at face value, these calculations would rule out the SD as one of the main channels 
to SN la because hydrogen has not been observed. Nevertheless, the issue is complicated 
because simulations predict that a large fraction of the contaminating hydrogen should be 
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moving at velocities below 1000 km.s~^, where it is mixed with Fe-peak elements that make 
its detection difficult because emission lines may blend with the many Fe and Co lines 
that appear during the nebular phase. On the other hand, observational constraints on 
the equivalent width of H lines in SN la spectra measured during the nebular phase put 
rough limits on the maximum amount of low-velocity hydrogen entrained in the SN ejecta. 



Mattila, Lundqvist & SoUerman (2005) give an upper limit of 0.03 M0 for SN 2001el, while 



Leonard (2007) quotes upper limits of 0.01 M© for SN 2005am and SN2005cf. 



The unambiguous observational detection of the companion star of the white dwarf 
(item b) above) is the most direct and effective way to constrain the progenitor scenario. 



A detailed search for such a star in the Tycho SNR was conducted by Ruiz-Lapuente et al. 



(2004), who claimed that a main sequence star with peculiar proper motion fulfilled all the 



criteria to be the companion. Subsequent studies of the properties of this main sequence 



star have not been so conclusive (Gonzalez- Hernandez, et al. 2009 Kenzendorf, et al. 2009). 



The imprint of the collision between the secondary star and the SN ejecta (item c) 
above) has also been analyzed in several works. In particular, the hole carved in the 
otherwise spherical ejecta by the shielding effect of the secondary is a source of asymmetry 



which could contribute to the observed diversity in SN spectra. Kasen et al. (2004) 



suggested a possible change of the peak magnitude with viewing angle of ~ 0.2 mag in B, 
comparable with the intrinsic dispersion of SN la light curves in this band. It was also found 
by the same authors that the hole is a source of polarization in the observed spectrum. Even 
though polarization is in general low in SN la spectra, it has been unambiguously detected 



and measured in surveys (see Wang & Wheeler 2008). Another observational signature 
of the collision could appear hours or days after the explosion as a persistent optical/UV 
emission from viewing angles looking down upon the shocked region, as suggested recently 



by Kasen (2010). Nevertheless, a recent search for this emission, as reported by Hayden et 
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al. (2010), did not give a positive result. All these effects would be detectable during the 



SN phase. There are comparatively few studies on the late stages of the evolution, once the 
ejecta start to interact with the ambient medium and the SNR phase begins. In a recent 



paper, Vigh et al. (2011) have analyzed the asymmetries introduced by holes with varied 



apertures in the geometrical properties of Type la SNRs using multi-D hydrodynamics. 
They suggest that the small asymmetries observed in the radius of the Tycho SNR can 
be interpreted as the imprint of the conical hole carved in the ejecta at the moment of 
the collision. In order to explain these asymmetries, the authors had to assume a very 



large angular width for the cone, ~ 90". According to Marietta, Burrows & Fryxell (2000) 



such large aperture is only possible if the secondary is a large red giant star, which would 
probably lose most of its weakly bound, H-rich envelope during the collision with the ejecta. 
Such a large amount of stripped hydrogen ~ 0.4 Mq should be visible in the spectra near 
maximum light and in the nebular phase, but this has not been reported in observations of 



normal Type la SNe. In a recent paper Lu et al. (2011) have invoked the SD progenitor 



model for Type la supernova as the cause for the prominent X-ray arc clearly visible in the 
SW quadrant of the projected disc of Tycho SNR. The increase of the observed brightness 



at the X-ray arc zone was interpreted by Lu et al. (2011) as the hallmark of the interaction 



between the supernova ejecta and the material stripped from the companion star shortly 
after the SN explosion. All this points to the SD scenario, in which the companion is an 
unevolved star, as a probable route to explain Type la SN explosions. 

The scarcity of studies of SNR evolution including the hole carved by secondary star 
at early times is the motivation for this paper. Past studies have highlighted the great 
interest to carry out simulations to elucidate if the hole will remain open or not after 



several hundred years of evolution (Marietta, Burrows & Fryxell 2000; Pakmor et al. 2008). 



Our main goal is to follow the hydrodynamic evolution of the hole region as the remnant 
expands into an homogeneous ambient medium (AM), and to infer possible observational 
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consequences. 

The calculations presented below do not refer to a particular Type la SN remnant, 
rather they attempt to outline the gross features that the presence of a companion star 
imprints to the long term evolution of the SNR. The secondary was assumed to be a main 
sequence solar-like star, a possibility that has been favored over red giant secondaries (see 
references in the discussion above). The choice of a companion mass of IMq was also 
motivated because that is the case for which the detailed study of [Marietta, Burrows fc ~ 



Fryxell (2000) provides more information (case termed HCV in their paper). In particular 
we want to check our estimation of the hole aperture and the velocity profile of the stripped 
matter inside the hole with that of Marietta and collaborators (described in their Sect. 4. 3). 
For example, they found that about half of the stripped material was confined in a cone 
within ^ 43*^ from the symmetry axis which can be taken as a rough measure of the size of 
the hole. Using AxisSPH we have obtained that ~ 40% of the stripped mass is within an 
angle of 40°. The slightly lower fraction of stripped mass within 6 = 40*^ from our simulation 
may be due to the differences in the assumed explosion energy of the supernova ejecta 
(~ 50% higher in Marietta et al.) as well as to the different nature of the hydrodynamic 
schemes used to carry out the calculations. Given the complexity of the calculations we 
have chosen to work with a unique scenario and to follow the interaction between the 
supernova ejecta and the companion star and later the surrounding AM. In this context, 
the choice of a solar-like star (a representative main sequence star) that is filling its Roche 
Lobe at the moment of the explosion is not unreasonable. 

Simulations were carried out using an axisymmetric SPH code developed recently. 



Garcfa-Senz et al. (2008). Because it was not possible to simulate the evolution over the 
large dynamic range in time between the collision and the fully developed SNR phase 
(minutes to thousands of years), we divided the process in three stages. First, we studied 
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the collision of the supernova ejecta with the secondary, spanning several hours. After that 
time, the interaction stops and both the ejecta and stripped material are in homologous 
expansion. At this point the secondary was removed from the calculation and the ejecta 
was stretched to a size of ~ 0.22 pc, roughly the radius of the system 28 yr after the SN 
explosion. This radius is large enough for the ejecta to swept an appreciable amount of AM 
material. Finally, a large region of uniform AM was introduced around the SN ejecta, the 
evolution of the SNR was followed until t~ 1000 yr, and the results were analyzed. 

This paper is organized as follows. In Section 2, we briefly describe the main features 
of the hydrodynamic method we use and describe the initial setting and the astrophysical 
scenario. Section 3 is devoted to the interaction between the supernova blast wave and the 
main sequence star, and to a comparison between our results and those of other authors. 
In Section 4.1 we describe the evolution of a spherically symmetric SNR as it propagates 
into the AM, as a benchmark to evaluate the asymmetric models that include the hole. In 
Sections 4.2 and 4.3 we study the evolution of these asymmetric SNRs. We conclude with a 
final discussion and a summary of our results. 



2. Hydrodynamic method and initial setting 

We carried out the simulations with the axisymmetric SPH hydrocode AxisSPH, 



described in Garci'a-Senz et al. (2008), which incorporates a new algorithm to solve 
the contribution of gravity in the axisymmetric SPH paradigm. The scheme calculates 
gravity by computing the direct interaction between any pair of particles, each of them 
approximated as a toroidal distribution of mass. To optimize computing time, the 
contribution of gravity was calculated only for the particles belonging to the secondary 
star. Thus, the problem was modelled as a free supersonic fluid, -the SN ejecta-, impacting 
on and passing through a gravitationally bound body, -the IMq companion star-. This 
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approach precludes us from estimating the degree of contamination of the secondary by the 
supernova material, but it is a reasonable approximation to model the structure of high 
velocity ejecta in homologous expansion. During the collision phase, the SN ejecta was 
represented by 10^ particles and the secondary star was simulated using 2 10^ particles. The 
initial model for both components was obtained by distributing the mass particles uniformly 
in a rectangular lattice and assigning them the mass required to reproduce the spherically 
symmetric density profile of SN ejecta in the homologous phase and a polytropic solar-like 
star of IMq. The spherical model used for the supernova ejecta was the deflagration model 



by Bravo et al. (|1996|) with kinetic energy at infinity of 8 10^° ergs. 



The equation of state (EOS) used to study the interacion of the supernova ejecta with 
the companion star was that of an ideal gas plus radiation. Shock waves were handled 
by adding an artificial viscosity term to the momentum and energy equations (see for 
example Rosswog 2009 for a recent review of SPH). In SPH, the linear and quadratic terms 
of artificial viscosity are controlled by parameters a and /3 respectively which we set to 
a = 1.5 and (3 = 3, slightly higher than the standard values to better handle with the 
hypersonic collision. This particular setting helps to prevent artificial particle penetration 
between the SN ejecta and the AM during the first stages of the collision. The same 



values for a and /3 were chosen by Herant & Benz (1992) to model the post-explosion 



hydrodynamics of SN 1987A with an axisymmetric SPH code, probably for the same 
reason. The center of the sun-like star was placed 2, 83Rq from the center of the explosion, 
at the point where the Roche-Lobe radius equals the radius of the secondary star. With 
that choice of parameters around 3% of the solid angle of the ejecta was intercepted by 
the secondary. This roughly corresponds to the HCV scenario considered by Marietta and 
collaborators (2000). 

Once the collision has come to an end, the companion star was removed from 
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the computational box and the ejecta structure was expanded homologously to a 
size of ~ 0.22 pc. The ejecta was then surrounded by an uniform AM with density 
Pam = 3 10~^^ g.cm~^ and size 4 pc x 8 pc. A large number of particles, Nam = 500, 000, 
evenly spread in a rectangular lattice, were necessary to encompass this volume. The 
structure of the ejecta was mapped with smaller number of particles N^j = 8800 in the 
innermost area, 0.22 x 0.44 pc, of the computational domain. As gravity plays a negligible 
role in this phase the algorithm used to calculate gravity was turned off. The EOS was 
changed to that of a perfect gas with 7 = 5/3. Once this basic spherically symmetric initial 
model was built it was adequately modified to host a hole with the angular size estimated 
from simulations carried out with AxisSPH as described in the next section. 



3. Description of the interaction between the supernova ejecta and the 

secondary star 



In the absence of direct observations, the effects of the impact of the supernova ejecta 
on the companion star must be addressed using numerical simulations. As stated above, 
the most complete numerical study to date is that of Marietta et al. (2000), who used 
an axisymmetric hydrocode to follow the dynamical phases of the collision process and 
estimated that ^ 0.1 — 0.2 M0 were lost by a solar-like companion. The amount of stripped 
mass from the companion star depends sensitively on the solid angle subtended by the 
secondary and to a lesser extent on the energy released by the explosion. The general 



picture of the process obtained by Marietta et al. was later confirmed by Serichol (2005) 



and Pakmor et al. (2008) using a three-dimensional SPH code. In the calculations of 



Serichol (2005), the orbital movement of the secondary was taken into account while it was 



not included in Pakmor et al. (2008) who conducted a resolution study of the event. The 



incidence of the orbital dynamics in the outcome of the collision is of minor importance for 



- 11 - 



magnitudes like the mass stripped from the companion, the size of the hole carved in the SN 
ejecta or the radial kick imparted to the companion star. However, other magnitudes are 
more sensitive to the orbital dynamics, especially those related to the mixing of the stripped 
hydrogen and the chemically stratified debris of the ejecta. In this case, the inclusion of 
the orbital velocity is crucial, because efficient mixing requires that the different materials 
coincide not only in physical space but also in velocity space. A quantitative comparison of 
several magnitudes of interest for the present work is given in Table 1. The differences in 
the stripped mass, the size of the hole, and the velocity kick between the well resolved 2D 
model A and the low resolution 3D models B and C of Table 1 are less than 20%, similar 



to the differences found by Pakmor et al. (2008) in their resolution study. Some of the 



discrepancies between the models shown in Table 1 come from the different prescriptions 
used to model gravity in two and three dimensions. 

The overall picture of the collision process can be seen in Figure 1, which shows 
density maps at different representative times for model A in Table 1. The last snapshot 
corresponds to t~ 5 h after the beginning of the collision. At this point, the interaction has 
ceased and the mass stripped from the secondary is O.IM©, in good agreement with models 
B and C, both calculated in 3D . The detailed temporal evolution of the amount of stripped 
mass in the three models from Table 1 is shown in Figure 2. About half of the stripped 
mass is removed from the secondary in an initial violent episode lasting around 250 s, and 
the remaining mass is removed more gradually during the interaction process. These two 
stages are usually referred to as the stripping and ablation phases [Wheeler, Lecar fc McKee 



(1975). Notice the brief episode of recapture of material which takes place around t= 210 



s in the three models depicted in Figure 2. It is not clear whether this feature, which is 



also present in the calculations of Pakmor et al. (2008), is real or a numerical artifact due 



to the simplicity of the criteria used to decide when a mass particle streams out from the 
secondary, namely that the actual velocity of the stripped material exceeds the local escape 



Fig. 1. — Density snapshots summarizing the evolution of the colhsion of model A in Table 
1. The first snapshot (top left) corresponds to t = 67s: the SN ejecta has compressed the 
frontal part of the envelope of the companion. At t = 190s (second snapshot, top right) 
the SN ejecta has already wrapped around the companion. The third and fourth snapshots 
(bottom left and right) correspond to times t = 340s and t=17922 s. The bow shock that 
is formed with the companion star at its apex can be clearly seen in all snapshots. The 
vertical dimension of the box changes from 2.2 10^^ cm (first sanapshot) to 2.5 10^'^ cm 
(fourth snapshot). 
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velocity. Such criteria neglect further hydrodynamic interactions of the stripped material 
being strictly valid only for material ablated just from the surface of the star. 

In Figure 3, we show the cumulative distribution of the SN ejecta and the stripped mass 
in model A as a function of 9, the angular distance to the axis defined by the centers of the 
SN and the secondary star (henceforth we will use 6 for this angular distance, and Oh = '2 
for the complete opening angle of the hole). In Figure 3 we can see how the hole in the ejecta 
geometry caused by the shadowing effect of the secondary is actually not empty, but filled 
with stripped material, basically H and He plus traces of heavier elements . Nevertheless, 
the mass of SN ejecta inside the hole is much lower than the mass stripped from the 
secondary. For angles smaller than 20° the region is almost devoid of SN debris, and it 
can still be called a hole. The slope of the mass distribution changes around 9 = 20° (see 
Fig 3), suggesting that the aperture of the cone carved by the secondary is ~ 40°, in good 
agreement to the calculations of Marietta et al. (2000) for a similar initial configuration. 
The amount of stripped gas within an angle 9h = 40° is ~ 0.03M -O.O4M0. 

In Figure 4 (left) we show the radial velocity profile for the particles initially belonging 
to the secondary at t~ 5 h after the explosion. All that stripped material is moving 
faster than the approximate escape velocity from the secondary, which is represented by 
the horizontal line at 520 km.s~^. As can be seen in Figure 4 (left) the stripped material 
is moving homologously, a feature which can be used to rescale the size of the system 
from hours after the explosion to years, when the interaction with the AM begins to 
affect the evolution of the remnant. In particular we have applied an homology ratio of 
5 X 10^ between the size of the ejecta at i ~ 5 h, the last calculated model of the interaction 
with the companion star, and its size at t ~ 30 yr, once ~ 2 10~^ Mq of AM have been 
swept by the SN ejecta. 

We have already seen that a little less than a half of the stripped material remains 
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inside the hole created by the shadowing effect of the secondary, while the rest is outside 
of this region, mixed with the SN ejecta. Past numerical simulations, [Marietta, Burrows &: 



Fryxell (2000), Pakmor et al. (2008), have shown that mixing between stripped material 



and SN ejecta only affects low-velocity Fe-peak elements, because the other elements 
synthesized in the SN are moving too fast to be mixed. The mass distribution in velocity 
space of the stripped material inside the hole is shown on the right panel of Figure 4. 
Most of the stripped mass is moving with velocities below 1000 km.s~^, with a peak in the 
distribution slightly higher than the escape velocity. This is also in agreement with the 



results of Marietta, Burrows & Fryxell (2000), who estimated a velocity of 823 km.s~^ at 
the half-mass point of the stripped material for a similar model. The profile of mass 
distribution in velocity space depicted in Figure 4 has a characteristic high and low velocity 
tails connected by a sharp line with a peak around log f (10^ km.s~^) = —1.3. While the 
high velocity tail is probably well resolved by the simulations the low-velocity region is not 
so well represented due to the difficulties to set an exact criteria to decide when a particle 
moving close to the escape velocity becomes unbound (for example it could lose velocity 
after colliding and be recaptured). The mass distribution in velocity space of the stripped 
material given in Figure 4 will be used in section 4 to set the initial model for the SNR 
evolution. 



4. Interaction with the ambient medium 



4.1. SNR models with spherical ejecta 



The large scale structure of SNRs is determined by the interaction between the SN 
ejecta and the surrounding ambient medium (AM). In this work, we consider a uniform AM 



with a density pam = 3 10 g.cm That choice is supported by Badenes et al. (2006) 



and Badenes et al. (2008) who have shown that the fundamental dynamical and spectral 



- 15 - 



Table 1. Main features of simulated Type la supernova ejecta and IMq star 



Model M, 



Mp 







» ej 

km.s"^ 



M. 



AM 



orb 



Mp) km.s ^ km.s ^ dee 



A (2D) 
B (3D) 
C (3D) 



1 2,83 
1 2,83 
1 2,83 



7500 0,04 
7500 0,04 
7500 0,04 



4,11 0,105 
4,11 0,09 
4,11 0,11 



232 



86 40 
72 43 
96 42 



Note. — When the orbit is included in the 3D calculation the width of the hole is slightly 



larger in the orbital plane ^// = 43 than in the ortogonal plane il± = 40 , Serichol (2005). 
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T3 

CL 

55 




Model A 
Model B 
Model C 



0.5 



2.5 



3.5 



1 1.5 2 

time (h) 

Fig. 2. — Evolution of the stripped mass for models A, B and C of Table 1. 
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Fig. 3. — Integrated distribution function of mass as a function of the azimuthal angle with 

the symmetry axis at t= 2840 s for model A. The red solid line corresponds to the angular 

distribution of accumulated stripped mass. Long dashed line in green is for SN material while 

short dashed line in blue is for the whole sample of supernova material plus stripped mass. 

For angles smaller than 20° the region contains only material removed from the companion 

star. 
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Fig. 4. — Left: radial velocity distribution of the material stripped from the companion in 
a cone of 6' = 20° at t= 17900 s for model A. The short horizontal line at 520 km.s^^ is 
the approximate escape velocity for the surface of an unperturbed 1 Mq polytrope. Right: 
Mass distribution of the stripped material in velocity space. 
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properties of Type la SNRs can be reproduced by models that assume a uniform AM. 



Also Badenes et al. (2007) showed that other kinds of AM (in particular, large wind-blown 
bubbles excavated by the accretion winds predicted by some Type la SN progenitor models) 
are inconsistent with the sizes of historical Type la SNRs. The adopted value for the density 



is representative for the warm phase of the ISM in our galaxy, Ferriere (2001), and close to 



the best fit value for Tycho. The typical size of the SNR at the time that a significant mass 

SMam 

MsN 



fraction / = ^rf^ has been swept by the blast wave is: 



R 



3fMsN 



3 



which gives radii of ~ 0.5 pc, 1.08 pc and 2.3 pc for f=0.01, 0.1, and 1, respectively. We 
have started the simulations of the SNR stage when the size of the SNR is Rq = 0.22 pc, 
corresponding to ^ 28 years after the SN explosion, to ensure that the swept-up mass is still 
very small. Such initial radii roughly corresponds to the size of the SN ejecta at the time 
when its outer density just equals the assumed density of the AM, pam = 3 10^^'^ g.cm~^. It 
corresponds to a value of f^ 0.1% in Equation (1). The starting time, t= 28 yr is, however. 



larger than that considered in the 2D calculation of Dwarkadas (2000) who made use of an 
expanding moving grid to increase the resolution during the very first stages of the SNR. 
As far as the initial swept-up mass is small enough the evolution of the remnant does not 
depend too much on the exact location of the contact discontinuity at the initial simulation 
time, 



Dwarkadas & Chevalier (1998). Taking / = 1 in Equation (1) we can obtain the size 
of the SNR when it enters the Sedov stage: R~ 2 pc. At that moment the age of the SNR 
would be a few hundred years. 

In order to benchmark the results obtained with the ejecta models modified by the 
collision with the secondary, we have also simulated the evolution of a SNR with spherical 
ejecta, quoted as model D in Table 2. We expanded the SN ejecta from its initial size to 
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a SNR age of t= 28 yr (outer radius ~ 0.22 pc). The density and velocity profiles after 
this homologous expansion are shown in Figure 5. At this age, the ejecta only fills a small 
fraction of the computational space, which, even in 2D, is a challenge for the numerical 
simulation. To build the initial model we have followed the method described in IVelarde ~l 



et al. (2006). The SN ejecta was represented by a small subset of N^tv = 8114 particles, 
evenly spread in a rectangular grid of (0.22 x 0.44) pc^, while the AM was represented 
with 'Nam = 562, 086 particles spread in a regular lattice of (4 x 8) pc^. The mass of the 
ejecta particles was adjusted to reproduce the density profile given in Figure 5. Although 
the total number of particles in the ejecta is small, the resolution is sufficient to follow the 
main features of the SNR evolution, including the development of the supersonic forward 



and reverse shocks (Velarde et al. 2006). 



The use of particles with different mass to reproduce the initial density profile warrants 
clarification because it could be the source of numerical artifacts which may distort the 
simulation. According to Figure 5 it is evident that the mass of the particles must change 
in a factor of the order of 10'^ in the neighborhood of the contact surface. Such large mass 
ratio may be a source of numerical troubles when the supernova particles and those of the 
AM mix. However, the real situation is not as bad as it may seem owing to the peculiarities 
of the mass of the particles in axisymmetric SPH codes. In the axisymmetric geometry the 

Table 2. Features of SNR models at t'=0.11 



Model Geometry Composition {6 < 20^ ) pam (g.cm ^) 



D Spheric Type la-SN 3 x 10" 

E Empty hole ■ ■ ■ 3 x 10" 

F Filled hole solar 3 x 10~ 
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local average of a physical magnitude A is calculated as < A >i= J2j 2m^^i^*i (where 
Vj is the distance of the particle to the Z-axis) from which < VA >j is conveniently 



estimated, Garcia-Senz et al. (2008). Therefore particles carry an effective mass m' = ^ to 
weight the interpolations. When the supernova material expands such effective mass is 
reduced. After one hundred years, when material begins to be mixed with the AM, the 
outer edge of the ejecta has already changed its radius by a factor ^ 5. When the ejecta 
reaches the limits of the system, that factor has grown to ~ 20. Due to the RT instability 
the supernova matter floats inside the AM whereas the less dense AM material falls down 
so that the initial contrast in effective masses is progressively reduced as the SNR evolves. 

Probably main effect of working with particles of different masses is that they increase 
the numerical noise which in turn serve as a seed of the RT instability during the first stages 



of the evolution. A similar conclusion was reached by Herant & Benz ( 1992 ) and Herant 



& Woosley (1994) who used an axisymmetric SPH code to simulate the post-explosion 
hydrodynamics of SN 1987A. The density snapshots depicting the evolution of the RT 
instability in these works are asthonishingly similar to these of ours, especially in the 
Herant & Woosley simulations (for instance compare their Figure 3 with our Figure 6 
below). What is reassuring is that they did also use particles with different mass but a 
quite different geometry for the initial lattice. 

To check more quantitatively the impact of using particles with different masses we 
ran a calculation in which the number of particles belonging to the ejecta was doubled 
(and their mass consequently halved). A second simulation was launched by doubling the 
number of particles everywhere so that resolution was a factor a/2 better. We have no 
detected significant differences in the evolution of the SNR. Apparently the main effect of 
doubling the number of particles in the SN region is to slightly delay the time at which 
the RT instability appears. Thus the impact of using particles with variable mass in the 
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simulations seems to be limited. It is a source of numerical noise that may distort the 
growth of the Kelvin-Helmholtz instability during the process of mixing, but we do not 
think it is seriously affecting the large scale features of the RT instability. It could also 
contribute to imprint the more filamentary look on the RT mushroom caps seen in the 
simulations presented below. 

The deceleration of the ejecta caused by the AM is equivalent to a local gravitational 
field pointing toward the center of the explosion, favoring the growth of the Rayleigh- Taylor 



(RT) instability in the dense layer between the forward and the reverse shocks. Chevalier, 



Blondin & Emmering (1992). In our simulations, the growth of the RT instability was 
induced by numerical noise in the initial distribution of the particles. As we will see, the 
development of the RT instability has a large impact on the long term evolution of the hole. 
Several snapshots of the density evolution of the SNR for model D are shown in Figure 



6 where times were given in normalized units t' = t/T being T defined as in Dwarkadas 



(2000): 



5/6 / N -1/3 



which for our parameter choice becomes T=250 yr. In the first snapshot, t' = 0.43 after 
the supernova explosion (t'sjm=0.32, henceforth we refer as t' and Vsim the normalized 
elapsed times since the explosion and from the beginning of the simulation respectively), 
the shocked ejecta still retain spherical symmetry. In the second snapshot, t'= 0.87, the RT 
instability has grown significantly, and the shocked ejecta begins to fill the interaction region 
between the reverse and forward shocks. The instability develops and gains complexity in 
the third and fourth images . In the last picture, the forward shock has already reached 
the limits of the computational domain, and the reverse shock, although still far from the 
explosion site, has gone back through more than the 95% of the supernova material. The 
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large plumes that develop at late times close to the symmetry axis are numerical artifacts 
due to the imposed axial geometry. The mass inside these plumes is very small, and it does 
not affect the outcome of the simulations in a significant way. The left panel in Figure 7 
summarizes the evolution of the spherical SNR at the moment when the blast wave have 
reached the limits of the computational domain. 

The onset of the self-similar Sedov stage depends on the AM density, kinetic energy of 
the ejecta and also on the precise profile adopted for the ejecta, [Dwarkadas fc Chevalier 



(1998). For the adopted values of these magnitudes in model D the swept-up mass becomes 



comparable to the ejecta mass at t'^ 0.92. Two dimensional simulations of SNR were 



(1992 


) and 


Dwarkadas 


(2000 



Dwarkadas (2000) investigated in some detail the interaction between the supernova ejecta 
and the ambient medium, either assuming constant AM as well as a CSM whose density 
decreases as r~^. These calculations were carried out using the finite-difference code VH-1 an 
assuming an exponential density profile for the ejecta. As our calculations were done using a 
SPH code with a different ejecta profile (see left panel in Figure 5) and a larger pam than in 



Dwarkadas (2000) a brief description of the evolution is warranted. During the self-similar 
evolutionary stage the average contact discontinuity radius evolves Rcd oc where the 
mean value for the exponent deduced from our simulations, s = 0.42 (see left panel in 
Figure 7) agrees with the theoretical value s = 0.4 expected for the Sedov solution. In the 



work of Dwarkadas (2000) the expansion parameter approach the analytical value after a 



time Vc^ 3 — 4. The normalized times used to find the expansion parameter in Figure 7 



go from t' = 0.9 to t' = 3.6 hence our results roughly agree to that of Dwarkadas (2000) 
but convergence to the theoretical value is a little faster. Such small discrepancy may arise 
from the differences in the ejecta profiles used in both simulations. There also differences in 
the morphology of the RT fingers, more filamentary in the SPH calculation, which extends 
to a larger radius in the unstable layer. On the other hand the initial phase of growth is 
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is slower in the SPH calculation owing to the lower resolution (in the VH-1 simulation a 
moving grid with a constant number of computational cells was used) and to the damping 
introduced by the artificial viscosity. 

The trajectory of the contact discontinuity is shown in the right panel of Figure 7, 
and it can be roughly fitted by a parabola, Rcd = 0.5 aco t'"^ assuming a deceleration 
a-CD = 0.2 pc in adimensional time units. . The growth rate F of the RT instability during 
the linear regime can be estimated using: 



A n acD , . 

where At is the Atwood number across the interface and n is the wave number . Once the 
swept-up mass has become comparable to the ejecta mass, at t' ^ 0.92 and Rcd — 2 pc, 
we get F = 0.22y^ for At = 0.5. After At' = 1, the e-folding growth factor, Cgf of a small 
perturbation is Cgf ~ 0.22^/n. Therefore, only perturbations with high wave number are 
in the non-linear regime at this time. However, in SPH small perturbations are usually 
damped by the artificial viscosity, so we expect that the RT instability will remain in the 
linear regime except at late times when the forward shock approaches the limits of the 
system. 

The location of the contact discontinuity follows the Rcoit') average trajectory, 
modulated by the effect of the RT instability. In Figure 8, we show the location of the 
contact discontinuity as a function of 9. To obtain the location of the CD, we recorded 
the position of the ejecta particle with the largest radius at each azimuthal angle. The 
fundamental features of this profile can be compared to those of the profile measured by 



Warren et al. (2005) for the CD in Tycho's remnant. The comparison, however, is not 
completely straightforward, because the profile observed by Warren et al. results from 
projecting the CD surface onto the plane of the sky, while our simulations provide a section 
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of this surface, relative to the explosion center, without projection effects. Another caveat 
stems from the axisymmetric hypothesis, which constrains the growth of the RT structures 



to have axial, ring-like, symmetry. In this respect Blondin & Ellison (2001) did not find 



large differences between two and three dimensional simulations for an EOS with 7 = 5/3. 
There is an increased amount of small scale structure in the 3D calculation as well as a 
slightly larger penetration of the RT fingers but the amount of turbulent energy in the 
unstable layer was similar. It is also worth to comment that for 7 < 5/3 (as it would be the 
case if there is particle acceleration at the shock front) the differences between 2D and 3D 



cases are larger, Blondin & Ellison (2001). 



The power spectrum (PWS) of the azimuthal distribution of the radial amplitude 
fluctuations of the contact discontinuity depicted in Figure 8 is shown in the leftmost panel 
of Figure 9. In the following section, we will compare these PWS to those resulting from 
ejecta models that include a hole. Both PWS shown in Figure 9 correspond to the last 
snapshot in Figure 6, but are calculated using only the upper {6 < 90°) and the bottom 
quadrants (90° < < 180°) of the computational space respectively. The PWS that we 



obtain is much flatter than in [Warren et al. (2005) (see their Figure 6), and it shows larger 
fluctuations. The mode with the largest power corresponds to n=l (or a wavelength in the 
azimuthal direction of 90°), which probably comes from the imposed symmetry between 
quadrants at t'sim=0. As expected, the PWS in different quadrants of the spherical ejecta 
simulation are very similar to each other. A linear fit to both spectra gives oc fc" with 
n ~ —0.5 , much lower than n ~ —1.5 obtained by Warren et al. (2005) for Tycho, although 



we remind the reader that this comparison is not straightforward. 
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Fig. 5. — Spherically symmetric density and velocity profiles of SN la ejecta at the time the 
simulation of the interaction with the AM starts. 
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Fig. 6. — Snapshots summarizing the evolution of the density during the Sedov- Taylor self- 
similar phase of the SNR corresponding to model D in Table 2. Adimensional elapsed times 
from the supernova explosion are t'=0.43, t'=0.87, t'=2.2 and t'=3.6. 
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Fig. 7. — Main features of the evolution of model D (left): the swept-up mass of the AM, the 
average radius of the contact discontinuity and the radius of the forward shock. A parabolic 
fitting for the locus of the CD is provided in the rightmost figure. The deceleration of the 
contact surface is a.cD — 0.2 pc in adimensional time units. 
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Radius of the CD as a function of the azimuthal angle 6 for models D, E and F of 



Table 2 at adimensional times t'=3.6, 3.9 and 3.9 respectively. 
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Fig. 9. — Power spectra of the CD locus shown in Figure 8 for models D (left), E (center) 
and F (right). Continuum lines are for the upper quadrant hosting the hole whereas dashed 
lines are for the lower quadrant. Fitting straight lines for the three models are also depicted. 
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Fig. 10. — Snapshots summarizing the evolution of the density during the Sedov- Taylor 
self-similar phase of the SNR once a. 6 = 20° cone of material belonging to the ejecta was 
removed from the simulation, model E of Table 2. Adimensional times are t'=0.38, t=0.81, 
t'=1.95 and t'=3.4 from the beginning of the simulation. 
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Fig. 12. — Intrusion of the supernova material toward the symmetry axis impelled by the 
RT instabihty (model E of Table 2). 
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Fig. 13. — Sketch of the coordinate transformation between the original cyhndric coordinates 
(r,z) to the one-dimensional projected coordinate (p) in a orthogonal direction to the line of 
sight defined by the viewing angle ip. 
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Fig. 14. — Projection of Model E in Table 2 in a orthogonal line to the line of sight for 
several viewing angles ip = 0°,20'^,45° and 60°. The normalized magnitude is depicted as 
a function of coordinate p in the projection line. 




Fig. 15. — Snapshots summarizing the evolution of the density during the Sedov- Taylor 
self-similar phase of the SNR corresponding to model F of Table 2. Adimensional times are 
t'=0.41, t'=0.86, t'=1.76 and t'=3.9. 
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Fig. 16. — Same as in Figure 12 but for model F in Table 2. 
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4.2. SNR from an initial ejecta with an empty hole 

We have first considered the situation in which the cone carved in the expanding 
supernova shell is initially devoid of any matter, model E in Table 2. This represents a 
limiting case for those models where the density in the hole would be much smaller than 
that of the AM at tsim=0 yr, being useful to compare to the spherical ejecta simulation 
described in the previous section and to the filled hole case described in the next section. 
Such hypothesis could adequately represent scenarios where the stripped matter allocated 
inside the hole is very small as it could be for instance the case of Tycho remnant where 



the total amount of stripped matter was estimated to be less than 8 10 Mq by Lu et al 



(2011). Several snapshots of the evolution of the remnant at different times are shown in 
Figure 10. Shortly after the interaction with the AM starts the hole begins to be replenished 
with AM material to finally be completely filled in a time At' ~ ^ , where Rq is the 
size of the ejecta at the beginning of the simulation. For model E in Table 2 that time is 
At' = 0.12. In the comoving frame moving at the average ejecta velocity the initially radial 
streamlines of AM diverts when they reach the edge of the hole at 6* ~ 20°. As the diffracted 
streamlines gain lateral momentum, they converge and compress in the neighborhood of 
the symmetry axis. Eventually, a steady state is reached in which the forward shock in this 
region lags behind the rest of the blast wave, distorting its spherical symmetry. Figure 11 
(upper rows) shows the velocity field in the upper quadrant viewed from the stationary 
center of our simulation space. At t'sim = 0.11, the hole is almost filled with AM material, 
which is compressed toward the hole axis. At t'sim=0.26, the velocity field in the hole has 
already aligned with the symmetry axis and the flow approaches a steady state. As time 
goes on, the RT instability develops, as shown in the second snapshot of Figure 10. At the 
outer edge of the hole the instability growth is particularly strong owing to the peculiarities 
of the hydrodyminamic flux in that region. Numerical and laboratory experiments with 



laser-produced plasma Kang et al. (2001) indicate that hydrodynamic instabilities and 
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vortex like structures form close to the symmetry axis when a supersonic flow of low density 
material goes through a stationary sphere made of higher density material. Moreover, in 
our case the presence of the hole also breaks the symmetry of the flow. As a result, the 
RT fingers at the edge of the hole grow stronger and project inside the hole volume at late 
times. This intrusion is illustrated in Figure 12, where we have overlayed a dashed line at 
— 20° for reference. Despite the strongly supersonic nature of the flow, this simulation 
suggests that the hole in the ejecta structure caused by the presence of the secondary star 
at early times could close or shrink appreciably over timescales of several hundred years 
due to the RT instability. 

In Figure 8 there is shown the location of the most distant particle of the ejecta a 
a function of the azimuthal angle. The profile can be compared to the that obtained in 
section 4.1 assuming spherical symmetry. As expected, the largest differences are found 
for 6 < 20°, but minor differences are also found at other regions in the upper quadrant 
where 9 < 90° (for example around 9 — 50°) while differences in the lower quadrant with 
90° < 9 < 180° are much smaller. This suggests that the initial asymmetry caused by the 
hole is affecting the hydrodynamic behavior of the remnant on a larger scale. In Figure 
9 (center), we show the PWS of the CD radius in the upper and lower quadrants of this 
model. At low wave numbers, the PWS in the upper quadrant is flatter than in the lower 
one and the slope of the linear flt is a little steeper. These results, although qualitative, 
suggest that the presence of a hole in the SN ejecta could be inferred by studying the 
geometrical properties of the CD in real SNRs. 

The dependence of the X-ray emissivity, (oc p^), of the shocked ejecta with the line 
of sight is quantitatively outhned in Figures 13 and 14. The procedure to make the 
projection of the remnant onto an orthogonal line to the line of sight is sketched in Figure 
13. For a given projecting line, with angle of sight V', the cylindric coordinates (r,z) 
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and density of the closest ejecta particle to the observer is recorded. The coordinates of 
these closests-to-observer particles with angle ip were projected into an orthogonal line 
to the line of sight (line p in Figure 13) and magnitude was reresented as a function 
of the unidimensional normalized variable p. The result of the projection for viewing 
angles ■0 = 0°, 20°, 45°, 60° is shown in Figure 14. For a viewing angle oi t/j — 0^ the 
observer is directly looking down into the hole. In this case the density profile around the 
origin is symmetric and there is a large density constrast between the hole region and its 
neighborhoods. As the viewing angle rises the density contrast goes down and the density 
gap moves to negative coordinates in the projection line. For viewing ang les ip < 60° the 
fingerprint of the 9h = 40° wide hole has already become very weak. Thus, on pure 
geometrical basis and taking ip ~ 45° as an (probably optimistic) upper limit to the viewing 
angle we would expect that ~ 25% {i/j = 45°) to ~ 10% {i/j = 20°) of Type la SNR could 
display inhomogeneities in the X-ray emission caused by the hidden hole. As the aperture 
of the hole primarly depends on the distance between the white dwarf and the secondary 
at the moment of the explosion, which is in turn related with the precise nature of the 
companion star, the detection of the hole could bring information about that point. 

4.3. SNR from an initial ejecta with a hole loaded with stripped material 

Several hours after the SN la explosion, the expanding ejecta is almost spherical except 
in a cone-like region with its apex located at the position of the secondary star at the 
moment of the explosion. This conical region is not empty, but loaded with H and He-rich 
material stripped from the companion. In order to model the effect of this material in the 
long term evolution of the SNR, we have included the basic features of the stripped gas in 
a 20° wedge around the symmetry axis. According to Figure 4, the stripped material in 
the hole region is moving homologously with a characteristic AM/Av profile which favors 
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the low velocity tail of the distribution. To incorporate the stripped mass, a 20° wide 
shce was removed from the SN ejecta. Then we took the same number of particles and 
assigned them solar composition and a velocity profile matching the one shown in the left 
panel of Figure 4. The mass of each particle was then modified to approximately follow 
the AM/A^; distribution given in the right panel of Figure 4, with the only constraint that 
the total amount of stripped material inside the hole equal the ~ O.OSSM© obtained in our 
detailed simulations of the stripping process. 

Several snapshots showing the evolution of the SNR are shown in Figure 15. On 
the whole, the SNR evolution looks very similar to the empty hole case explored in the 
previous section. Nevertheless, the addition of the low- velocity material stripped from the 
secondary introduces a few significant differences. For example, the flow inside the hole 
gets aligned faster with the symmetry axis, as can be seen by comparing the upper and 
lower rows in Figures 11. At the outer edge of the hole, there is a larger stirring effect of 
the ejecta material as it interacts with the low density but high velocity component of the 
stripped matter. As shown in Figure 15, this leads to a stronger development of the RT 
instability in that region. At t'=1.76 the RT finger around 9 ~ 25° is bigger than in the 
the empty hole case, and its outer edge is close to the forward shock. The high density 
but low velocity component of the stripped material does not have a major impact on the 
SNR geometry, partially due to the damping introduced by the artificial viscosity in our 
SPH simulations, which suppresses the growth of instabilities close to the center. To follow 
the hydrodynamical evolution of the innermost region of the hole with sufficient detail, a 
much higher resolution study should be conducted. As in the previous section, the strong 
development of the RT instability at the edge of the hole leads to the intrusion of some 
supernova ejecta into the hole (see Fig 16), but to a lesser extent. The final fate of the hole 
is not clear. Although the hole may never close completely, our simulations do not exclude 
its partial closure by hydrodynamic instabilities at late times. In any case, the simulations 
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presented here suggest that the hole will remain open during several centuries, distorting 
the geometry ofthe CD in historical Type la SNRs. 

The radius of the CD is at t'=3.9 is shown in Figure 8. The profile is similar to the 
empty hole case, with minor differences. As discussed above, the hole seems to have closed 
more in the ejecta profile without stripped mass than in the one including it. The radial 
amphtude of the CD at ^ = 25° is bigger owing to the larger vertical excursion of the RT 
mushroom in that region. These pecuharities are also present in the PWS of the angular 
distribution of radial amplitude of the CD, as shown on the right of Figure 9. The PWS 
in the upper quadrant is similar to that of the empty hole case but the features are more 
pronounced and the slope of the linear fit steeper. The PWS in the lower quadrant was 
similar to that of spherical and empty hole calculations, as expected. As in the precedent 
section we have projected the supernova ejecta onto a line orthogonal to the line of sight 
(Figures 13 and 17) to estimate limiting viewing angles to detect the hole. The results were 
similar to that of the empty hole case. 

5. Conclusions 

The imprint of the secondary star in the long term evolution of Type la SNRs has 
been studied by means of numerical hydrodynamic simulations. We simulated first the 
interaction between the spherically symmetric SN ejecta and the nearby companion star 
in the standard SD scenario for Type la SNe. We used the results of these simulations to 
set the initial conditions for the later phases of the SNR evolution. Our calculations were 
carried out using an axisymmetric SPH code, which makes it easy to keep track of the 
several components of the fiuid, namely the supernova material, the AM and the material 
stripped from the envelope of the companion star during the initial interaction. 
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The results of the first phase were in basic agreement with those of [Marietta, Burrows &: 



Fryxell (2000), Serichol (2005), and Pakmor et al. (2008) for the chosen initial configuration 



of the binary system at the moment of supernova explosion (summarized in Table 1). 
Of special relevance for our study were the angular amplitude of the hole opened in the 
supernova debris by the shielding effect of the companion, sun-like star, the amount and 
velocity profile of the stripped material and its distribution in velocity space in the hole 
region. 

Once the initial interaction ended, an homologous transformation was applied to set 
the initial conditions for the second phase and a large volume of AM was incorporated to 
the computational domain. The ratio between the volumes encompassed by the AM and 
the supernova debris at t=0 yr was ~ 6000. Such a huge value highlights the difficulty 
to carry out a full three dimensional study of the phenomena with sufficient resolution. 
The existence of a symmetry line makes it possible to handle the evolution of the remnant 
using an axisymmetric hidrocode so that the achieved resolution was enough to grasp the 
main features of the of the interaction of the SN ejecta with their surroundings. Three 
calculations were carried out, assuming complete spherical symmetry in the ejecta, an 
empty 9h = 40^ hole due to the presence of the companion, and a similar hole filled with 
stripped material from the secondary star. Table 2. 

Our simulations show that the hole carved in the ejecta affects the long term evolution 
of the SNR. We have seen that hydrodynamic instabilities at the edge of the hole trigger 
the intrusion of material toward the symmetry axis, especially in the empty hole model. 
When the hole is filled with material from the companion star, this intrusion was somewhat 
suppressed. We conclude that the conical hole will remain open during a long time, probably 
longer than the typical age of historical SNRs, but its closure over longer timescales is not 
ruled out. 
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We have characterized the influence of the hole on the geometrical properties of 
the contact discontinuity separating the supernova material from the AM. Within the 
limitations of our study, we showed that the CD has different geometrical properties in 
the hemisphere that contains the hole. The hole seems to influence the dynamics of the 
SNR over angular scales that are larger than its size. The Fourier analysis of the angular 
distribution of the radial fluctuation of the CD suggests that the power spectrum of the 
CD surface should have a steeper exponent close to the hole than away from it, but this 
statement should be conflrmed by a full 3D study with sufficient resolution. 

Another point of interest is that hydrodynamic instabilities at the outer edge of the 
hole develop faster than the average growth rate in the RT unstable region of the SNR, and 
can come close to the forward shock. The extension of the RT instability layer has been 



addressed in several works, Blondin & Ellison (2001), Wang & Chevalier (2002) where two 



mechanisms were identified as agents which could significantly enhance the penetration of 
the RT fingers. One is the reduction of the effective adiabatic 7 value used in the EOS due to 



particle acceleration at the forward shock, Blondin & Ellison (2001). In this case the larger 



extent of the region subceptible to the RT instability would be more or less homogeneous, 
thus affecting the whole unstable shell. The second mechanism invokes the existence of 
isolated dense ejecta clumps in core-collapse SNR to explain the observed protusions in the 



Vela supernova remnant, Wang & Chevalier (2002) and would be much more localized than 



the first one. According to our calculations the existence of a void in the supernova ejecta, 
either filled or not with stripped matter from the former companion star, may also leads to 
a larger development of the RT fingers in the outer edge of the hole region. This suggests 
that, if observed, the existence of isolated RT structures with anomalous size in putative 
Type la SNR could be an indication supporting the SD escenario. 

Nevertheless the limited size of hole would make its detection difficult. According 
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to the geometrical analysis given in Sections 4.2 and 4.3 the chance to detect the hole is 
small for viewing angles (with respect to the symmetry axis) i/j > 30 — 40° and completely 
negligible for ip > 60''. For a random SNR distribution that means that only one among five 
or six SNR arising from Type la SN explosions which hosts a ~ 40° hole may be detected 
this way. Such probability is lower /higher for smaller/wider holes, which in turn depends 
on the particular features of the binary system where the explosion of the white dwarf took 
place. 

The novel calculations described in this manuscript confirm that multidimensional 
hydrodynamics can be safely applied to the study of the long term evolution of SNR after 
the explosion of a SN la in the single degenerate paradigma. Given that simulations were 
much time consuming we have focussed on a particular combination of parameters of the 
binary system hosting the explosion and used an uniform ambient medium to explore the 
consequences on the long term evolution of the remnant. Therefore our analysis does not 
apply to a particular SNR although we expect that the gross features of the phenomena 
were captured by our simulations. 

We thank the anonymous referee for the constructive critical comments that much 
improved the presentation of the manuscript. This work was funded by Spanish MICINN 
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